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Abstract 

We present a closed description of the charge carrier injection process from a conductor into an 
insulator. Common injection models are based on single electron descriptions, being problematic 
especially once the amount of charge-carriers injected is large. Accordingly, we developed a model, 
which incorporates space charge effects in the description of the injection process. The challenge 
of this task is the problem of self-consistency. The amount of charge-carriers injected per unit 
time strongly depends on the energy barrier emerging at the contact, while at the same time 
the electrostatic potential generated by the injected charge- carriers modifies the height of this 
injection barrier itself. In our model, self-consistency is obtained by assuming continuity of the 
electric displacement and the electrochemical potential all over the conductor /insulator system. 
The conductor and the insulator are properly taken into account by means of their respective 
density of state distributions. The electric field distributions are obtained in a closed analytical 
form and the resulting current-voltage characteristics show that the theory embraces injection- 
limited as well as bulk-limited charge-carrier transport. Analytical approximations of these limits 
are given, revealing physical mechanisms responsible for the particular current-voltage behavior. 
In addition, the model exhibits the crossover between the two limiting cases and determines the 
validity of respective approximations. The consequences resulting from our exactly solvable model 
are discussed on the basis of a simplified indium tin oxide/organic semiconductor system. 

PACS numbers: 73.40.-c, 73.40.Lq, 72.80.Le 



Electronic address: yugenen@tgm.tu-darmstadt.de 



1 



I. INTRODUCTION 



Once a conductor forms contact with an insulator, an energy barrier is formed between 
the two materials, which impedes the charge- carrier injection into the insulator. Although 
this injection barrier in general sways charge transport through the conductor /insulator 
system, only the two limiting cases of very low or very high injection barriers are often 
considered. 

For low injection barriers, one expects the contact to be ohmic, meaning that the contact 
is able to supply more charges per unit time than the bulk of the insulator can support, 
n this case, a space-charge region is formed and the electric field at the interface vanishes 
l|. Because excess charge-carriers dominate charge transport in insulators, one observes 
a space-charge-limited current (SCLC) density of the form j ~ V 2 /L 3 (in the absence of 
charge-carrier traps), where L is the sample thickness and V is the applied voltage. The 
current-voltage characteristic (IV-characteristic) is determined by the bulk properties of the 
material with no influence of the contact properties j^, [3, 0]. 

For high injection barriers, one anticipates the injection rate across the conduc- 
tor/insulator interface to dominate the IV-characteristic of the system. The models to 
describe injection are the Fowler-Nordheim (FN) tunneling model [5J or the Richardson- 
Schottky (RS) model 5|, |6| for thermionic injection. The FN model describes the charge- 
carrier injection by tunneling through a triangular barrier into an unbound continuum of 
states. The current takes then the temperature independent form: 



V 2 

j(V) = B j2 eX P 



3heV 



(1) 



where e is the elementary charge, m e ff being the effective mass of a charge-carrier, H is 
Planck's constant, B is a constant and A is the height of the barrier. The RS model on the 
other hand describes charge injection as a thermally activated hopping over the potential 
barrier, where barrier lowering due to the superposition of the external electrostatic potential 
and the image-charge potential is considered. The RS model predicts the IV-characteristics 
to follow: 

, W . C r*exp[-A] e xp[±^]. (2) 

Here, T is the temperature, k is the Boltzmann constant, C represents the effective Richard- 
son constant, e is the relative permittivity of the insulator and eo is the permittivity of 
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vacuum. 



An alternative description of the injection process is given by the dri: 



t-diffusion theory 
8fl. In this homo- 



involving electron-electron interaction in a mean-field approximation 0. 
geneous continuum model widely used for the description of conventional crystalline semi- 
conductors, the steady state current density exemplarily given for holes in one dimension 
is: 

j = efi s p s (x)J : ' s (x) - eDp' s (x) = const, (3) 

where fj, s is the hole mobility, p s is the hole density, T s is the electric field in the insulator, 
and D is the diflusivity of holes. T s and p s are coupled by the Poisson equation. The drift- 
diffusion equation in combination with the Poisson equation involves space-charge effects, 
but meets the problem of self-consistency in the boundary conditions. The electrostatic 
potential generated by the injected charge-carriers modifies the injection barrier, but on 
the other hand the amount of charge-carriers injected per unit time depends strongly on 
the barrier height. Therefore, a proper description of the injection process by means of the 
drift-diffusion theory has to involve both, the insulator and the conductor side of the system. 
This can be seen in analogy to a strongly asymmetric jm-junction. 

Different attempts to numerically describe the inj ection process in insulators with ac- 



c 0Unt of t fce space cllarg e were m ade re ce ntly fl fl U In transitlve FMmutons m 
one- and two-dimensional geometries by Christen et al. Ill] both injection-limited and space 
charge-limited regimes were presented. However, the injection description was still not self- 
consistent imposing some phenomenological values for the field and the carrier density at the 
contact. A comprehensive one-dimensional numerical model accounting for space-charge ef- 
fect was developed by Tutis et al. [9] which comprised the hopping transport and the tunnel 
injection from electrodes. Unfortunately, this sophisticated numerical tool does not allow 
analytical fitting of current-voltage characteristics which gives insight in major mechanisms 
controlling injection. In the work of Neumann et al. [10!], a self-consistent numerical treat- 
ment of the injection and transport processes in a conductor /insulator /conductor device was 
presented in which continuity of the electrochemical potential and the electric displacement 
was assumed everywhere in the system, in particular at the contacts. The conductor and the 
insulator were characterized by their specific density of state (DOS) distributions. In present 
work we obtain an exact analytical solution for the injection across the conductor /insulator 
interface in a self-consistent manner. To focus on the consequences of the self-consistent 



treatment, we discuss our results on the basis of a simplified indium tin oxide/organic semi- 
conductor system used as hole-injecting contact in organic optoelectronic devices, shifting 
the comparison with experimental data to a future time. 



II. THE MODEL 

Let us consider a single conductor/insulator interface located at a position x = 0. The 
conductor is supposed to extend over the semi-space with x < 0, whereas the insulator covers 
the semi-space with x > 0. All energies are measured downwards with respect to the top 
of the valence band in the conductor, in order to account for hole transport in the system. 
In the following sections, the theoretical models describing the insulator and the conductor 
are introduced. 



A. Electrode 

The conductor electrode is characterized by its DOS as a function of the energy E. For 
our purpose it is sufficient to assume the free electron approximation, in which the DOS 
function reads: 



1 (2m eff 



*w = s? -iT) (4) 



where m e ff is the effective mass in the conductor. In the Thomas- Fermi approximation 
one can deduce the electrochemical potential k c of the conductor as a function of the spatial 
coordinate x, 

Kc(x) = -^—(3n 2 p c (x)) 2 / 3 + e0(x), (5) 
2m eff 

where p c {x) is the hole density in the electrode and 4>(x) is the electrostatic potential. In 
general, the electrochemical potential k(x) relates the steady state current density j with 
the charge-carrier density. For a one-dimensional geometry, the current remains constant 
across the whole space and j is given by the conductivity a and the derivative of k(x) 

= _adK(x)_ 

e dx 

(notice the direction of the energy axis). The conductivity of a conductor o c = t[i c Poo can 
be expressed in terms of the hole mobility fi c and the hole density p^ in the valence band 
at an infinite distance from the conductor /insulator interface. 



Since charge-carriers are transferred from the electrode to the insulator, a space-charge 
region emerges near the interface which modifies the electric field ^(x) in the conductor 
according to Gauss law, 

r c (x) = —6p(x), (7) 
e c e o 

where e c is the relative permittivity of the electrode. In Eq.flZj) 8p(x) = p c (x) — is the 
excess hole density. However, charge- carrier densities in degenerate conductors are rather 
high and consequently, the value for the excess hole density is small in comparison with the 
background hole density, |5p(x)| <C Poo- Hence, the well-known linearized Thomas- Fermi 
approximation [g] can be applied, leading with Eqs.([MZD to a differential equation for T c : 

l 2 TF K{x)-T c {x) = -±- (8) 

with 



V 3 e 2 Poo 

being the Thomas-Fermi screening length, defining the typical length scale of the system. 
Here is the electrochemical potential at an infinite distance from the conductor /insulator 
interface. 

Since space-charge zones in conductors are of finite thickness, gradients of T c {x) have to 
vanish at an infinite distance from the contact. Therefore, the solution for the electric field 
reads 



TAx) 



^c(O) - - 

Or 



e x/l TF + J_ ; (10) 



where the electric field in the conductor at the conductor /insulator interface, ^(0), is the 
only unknown quantity. 

The validity of the used approximations has to be reviewed critically. An exact quantum- 
mechanical theory of the inhomogeneous electron gas accounting for the ionic lattice of the 
material and electron correlations, gives a comparable scale for the electric field penetration 
in conductors and demonstrates the usability of the uniform positive background model 
for simple metals Q]- While in typical metals the Thomas- Fermi screening length 
is about lA making the application of the quasi-classical Thomas-Fermi approximation, 
Eq.flHJ), inappropriate, substantially larger Thomas- Fermi screening lengths can be found 
in degenerate semiconductors such as transparent conducting oxides used in optoelectronic 



devices 
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B. The insulator 



Similarly, the insulator is characterized by a DOS function, g s {E), describing extended 
states in which charge transport takes place. In contrast to the conductor the electrochemical 
potential, k s , is situated well above the top of the band of extended hole states. In other 
words, the insulator is supposed to be non-degenerate. Introducing a band-edge means that 
the DOS function g s (E) = when E < 0. Hence, the density of holes in extended states 
can be calculated using Boltzmann statistics: 

p.(x)= J 9s (E — A — Kpo) exp ( Ks(x) ~ e ^ x) ~ E ) dE, (11) 

— oo * ' 

where the energy scale has been adjusted to the top of the valence band in the conductor. 
The injection barrier A is defined as the energetic difference of the top of the extended state 
distribution, g s {E), to the electrochemical potential in the conductor at an infinite distance 
from the interface, k^. The electrochemical potential of the insulator is then given by 

k s (x) = kT\n + A + /too + e0(x), (12) 

where the quantity 

oo 

M = J g s (E)exp(-E/kT)dE (13) 

— oo 

can be understood as the effective total density of states available in the insulator at a given 
temperature. We note that the temperature dependence of M becomes weak in the case of 
a narrow-band insulator. 

Since band-gap energies are much larger than kT, thermal excitation of a charge carrier 
from the valence band to the conduction band of a typical insulator is virtually impossible. 
Charge- carriers contributing to the electrical current are therefore excess charge-carriers 
injected from the conductor and their total density has to appear in Gauss law, 

T' s {x) = —p s (x), (14) 

with e s the relative dielectric permittivity of the insulator. Eqs.(j5]). (|12p and (|14p with 
a s (x) = efi s p s (x) lead to a nonlinear differential equation for the electric field J r s {x), 

kT i 

—r s \x)-F s (x)r s (x) = — 3 — (is) 
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where /i s is the hole mobility in the insulator. The same result is obtained by employing the 
drift-diffusion model, Eq.Q, and the Einstein relation, relating /i s and D [6]. Introducing 



the following dimensionless quantities, 

X 



I 



-x. 



Itf 
elrF 
kT 



T 

J Si 



e 2 / 3 



TF 



[i s e s e (kT) 2 - 



Eq. (fl5l) converts into a dimensionless form, 

F' s '(X)-F' s (X)F s (X) + l = 0. 



(16) 
(17) 
(18) 

(19) 



The solution of Eq. (|T9|) must be separately formulated for the two cases of thermal 
equilibrium and an applied steady-state current. In equilibrium, the dimensionless current 
density i vanishes and Eq. (|T9|) can be integrated in elementary functions. The first integra- 
tion results in: 

F' S (X) - l -F*{X) = A (20) 

with an arbitrary constant A. At an infinite distance from the contact the field and its 
derivative (i.e. the charge carrier density) vanish, thus, A = and a solution for the electric 
field reads, 

Fs(X) = - ^ ( °L . (21) 



l-F s (0)X/2' 

In Eq. (1211) the field at the contact in the insulator, F s (0), is the only unknown quantity. 
Considering a net current density, a general solution of Eq. (I19p is known in terms of Airy 



functions Ai and Bi 
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19j: 



/3 1/3 Ai' [{l/2) 1/ \X + Ci)] + C 2 Bi' [(,/2)V 3 (X + C,) 
{ } 1 Ai[(6/2)V3(x + C 1 )]+C 2 Bi[(V2) 1 /3(X + C 1 )] 



(22) 



where primes denote derivatives of Airy functions with respect to their arguments, C% and 
C*2 are unknown constants. When X — > oo the gradient of the charge carrier density has to 
vanish so that, 

F S (X)F' S (X) = i. (23) 

Eq. (l23l) demonstrates that, in presence of a constant current, the magnitude of the electric 
field F S (X) rises asymptotically since the charge carrier density ~ Fg(X) vanishes. Con- 
sidering the asymptotic behavior of the solution of Eq. f|T9|) it is convenient to account 



explicitly for the sign of the current density i. Since we assume injection of holes from the 
left semi-space to the right one, i together with j must be positive. This imposes asymptotes 
of the field and the hole density 



F S (X) ~ VX, Ps ~ 1/VX, (24) 

respectively, resembling the characteristic behavior for SCLC [1]. Taking into account 
asymptotic properties of Airy functions 19] it is easy to establish that, to satisfy the asymp- 
totic conditions for the field, Eqs. (1231) and (|24l) . the constant Ci must equal zero. Then, 
the general solution, Eq. (1221) . transforms to 

Ai' \(t/2) 1 / 3 (X + C 1 ) 



F S (X) = -2(./2) 1 / 3 - 



(25) 



The constant C\ must be determined from the boundary conditions derived in the next 
section. 



C. Self-consistency and boundary conditions at the contact 

The equations of the electric field distributions in the conductor and in the insulator 
have to be solved self-consistently. Assuming steady-state, self-consistency is achieved by 
adjusting the integration constants (i.e. F s (0), F c (0) and C\) with respect to the Maxwell 
equations and continuity of the electrochemical potential: 

eJ-(x) = continuous, (26) 
k(x) = continuous. (27) 

Eq.( l26l) requires that the electric displacement in the insulator at the interface, e s F s (0), 
is equal to e c F c (0). This holds as long as no interface char ges or dip ole- layer exist at the 
contact, being problematic in many systems considered here [2l|, l22j|. Yet, for simplicity 



and without loss of generality Eq. fT26l) is assumed, leading with Eq. fTSTl) and Eq. fTlOl) to 
a nontrivial boundary condition relating the electric field and its first derivative in the 
respective media at the contact. In the insulator the following boundary condition holds: 

-^(°) - Pf - ln («^(°)) = "*~ > (28) 

Cc "<-/ jJ, c Poo 
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where the dimensionless constant a is defined as 

^ _ e s e kT _ 3 e s kT 

This boundary condition contains parameters of both media, being specified in the bulk of 
the respective material. 



III. PHYSICAL AND NUMERICAL ANALYSIS 

In this section the solution of the injection problem using the presented model is discussed. 
We distinguish between the equilibrium condition, where space charge zones are formed as a 
consequence of diffusive charge carrier transfer, and the case of a steady state current, where 
charge carriers are driven through the system by a time invariant external electric field. 

As an example for an insulator, an organic semiconductor can be considered. Organic 
semiconductors show many typical characteristics of insulators like relatively large band gaps 
up to 3 eV and hence, the absence of intrinsic charge carriers. However, it is well established 
that disordered organic semiconductors possess a Gaussian DOS |23j] compromising the 
applicability of Eq. ffTTl) - organic semiconductors are degenerate systems, the tail states 
acting as charge carrier traps. As a consequence, one has to distinguish between trap states 
and transport states and Fermi statistics has to be considered. Yet, for weak disorder, 
the Gaussian width is small and in the limiting case of a vanishing disorder, charge-carrier 
trapping in tail states is negligible and Boltzmann statistics is valid. 

In organic light emitting diodes (OLEDs) or field effect transistors (OFETs) organic 
semiconductors are contacted with metals like Au, Ca, Al or transparent conducting oxides 
like indium tin oxide (ITO) to allow for charge carrier injection in the otherwise charge- 
carrier free organic semiconductor. While in metals the Thomas-Fermi approximation is 
disputable due to the prevailed low screening length of lA, the characteristic scale Itf varies 
in ITO, depending on doping [jj], 16, 17], from 2AA to a few nm. Therefore, we choose in 
our self-consistent consideration an ITO electrode and assume that it can be described in 
lerms of the Thomas- Fermi approximation. ITO is typically employed as anode in OLEDs 
8|, since it provides a decent conductivity and a sufficient high workfunction (5eV) to allow 
for efficient hole injection while being transparent in the visible range of the optical spectrum 
to enable light outcoupling. 
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organic 


ITO 


N 


e s 


Us 


Poo 


m eff 




e c 




Itf 


(cm -3 ) 




(cm 2 /Vs) 


(cm -3 ) 


(m e ) 


(eV) 




(cm 2 /Vs) 


(A) 


10 21 


3 


10~ 4 


10 20 


0.35 


0.225 


9.3 


30 


8.6 



TABLE I: Typical material parameters for an organic semiconductor and ITO 
parameters are deduced assuming T = 300i"C. m e is the electron mass and a 
determined by parameters of both materials. 



fl fl Q 



5.6 • 10 -3 



f/. The 

is 



Due to the importance of charge carrier injection for the device performance of OLEDs 
and OFETs, the description of the involved charge carrier injection process has been ad- 
vanced in recent years considering surface recombination of charge carriers at the interface 
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2811 . Even more 



251 ] or stochastic hopping in a surface-barrier potential 
suitable for organic semiconductors microscopic models account for the mobility of electrons 
within and transfer between molecular strands and for the interaction of electrons with 



molecular vibrational modes 29(. However, all these models work within the single electron 
picture, so far it concerns the injection process, which means that interaction between in- 
jected electrons is not incorporated and, thus, space-charge effects on the injection are not 
properly taken into account. 

From now on it is assumed that the material specific quantities of the organic semicon- 
ductor and ITO adopt the typical values given in Table [B Thereby, the injection barrier A 
is given by the energetic difference between Koo and the band edge in the organic semicon- 
ductor and, thus, is determined by both media. Changing the value of the barrier height 
while leaving the electrode unchanged can therefore be understood as considering a different 
organic semiconductor. 

Specifying the material parameters for the insulator and the conductor leads to substan- 
tial consequences for the boundary condition given in Eq.( l28i) . l is multiplied by a small 
factor, a(/jL a / fjL c )(Af /Poo) ~ and hence, the boundary condition does not depend directly 
on i in most practical cases. Neglecting ta(/i s //i c )(JV r /Poo)> Eq. (l28l) can be reformulated to 



p s (0) = A/exp [(e s /e c )F s (0) - A/kT] 



(30) 



p s (0) being the density of injected holes at the contact. As a consequence, the dependence 
of F s (0) on the current is only due to Eq. ffl9|) . Apparently, there exists a barrier variation 
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oc F s (0) leading to the definition of an effective injection barrier A e g of 



A eff = A - kT(e s /e c )F s (0) 
= A-eF c (0)l TF . 



(31) 
(32) 



The modification of the injection barrier corresponds to the amount of energy a charge 
carrier gains (or loses) in the electric field at the electrode side of the interface. This change, 
as is seen from Eq. (1321) . may be positive or negative, since the space-charge region at the 
interface represents a potential barrier itself. 

A. Equilibrium 

Using the analytical solution (12 ip valid for thermal equilibrium, the boundary condition 
([2"g]) can be simplified to a transcendental relation for F s (0): 



F s (0) together with Eq. (l21l) yields the solution of the electric field distribution in thermal 
equilibrium. While the charge carrier mobilities of the respective material determine the time 
needed to reach equilibrium they apparently do not influence the final equilibrium electric 
field distribution. Eq. (f2"Tj) and Eq. (l3"3l are generally valid for non-degenerate systems. 

The solution for the dimensionless electric field is shown in Fig{TJ where barrier-free injec- 
tion (A = 0) from an ITO electrode in an organic semiconductor is considered. The solution 
can be interpreted as follows. Holes diffuse from the electrode into the organic semiconduc- 
tor. This results in a negative electric field causing a drift current opposite but equal in 
absolute value to the diffusion current, so that the net current is zero. As a consequence, 
a space-charge zone is established in the electrode and in the organic semiconductor, which 
in total is neutral. This space-charge region is very thin in the electrode but extends far 
into the organic semiconductor. The weak decay of the electric field (|F<,(A)| ~ 1/^) in the 
organic semiconductor is due to a missing charge- carrier background in the insulator. The 
distribution of the electric field changes when a non- vanishing injection barrier is introduced. 
Due to an impeded charge carrier injection, the electric field is solely reduced close to the 
interface, leaving the field far in the organic material invariant. 

Depending on the material parameters chosen, the parameter a, Eq. (1291) . may be very 
large or very small (here: a = 5.6 • 10~ 3 ). In these limits the nonlinear equation (133]) can be 




(33) 
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FIG. 1: Distribution of the electric field F in units of Fq 



kT 

el TF 



in equilibrium (j = 0) and barrier- 



free case (A = 0) as function of the coordinate X in units oflxF- The inset shows the discontinuity 
of the electric field at the interface. 



solved approximately. If a << 1, two characteristic regions arise depending on the relation 
between A and Ao = kT\n(2/a) (here: Ao = 0.15 eV). If the barrier height is so small that 
A << A , F a (Q) is given by: 



F,(0) = -(Ae c /e s ) 1 - (1 + A/2)" 1 ln(A e c /e s ) 



with A = ln(2/a) — A/kT » 1. If the barrier is so large that A >> A then 



Fs(0) = -\ -exp 
a 



2kT) 



(34) 



(35) 



The latter relation is also valid for a » 1 and arbitrary As. Eq.( !35l) demonstrates the 
exponential suppression of the electric field within the space-charge region by the injection 
barrier, since less charge carriers are transferred across the contact. 

The integration of the e.eetde field leads to a voltage V eqm to maintain e q uihb r iu,„. The 
equilibrium voltage V equi is often referred to as the contact potential [5J . Since the equilibrium 
field decreases with 1/X, the integration over the entire infinite semi-space diverges so that 
a cut-off length L has to be introduced. This results from the fact that a perfect insulator 
excluding any intrinsic charge carriers is considered. The analytical expression for V eqU i 
reads, 



equi 



kT 

e 



-F s (0)-21n 1 



F s (0)L 



21 



TF 



(36) 
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where F s (0) results from the boundary condition in equilibrium, Eq. (!33l . The existence 
of such an equilibrium voltage has practical consequences for organic electronic devices. In 
organic photovoltaic cells, the contact potential is known to reduce the open-circuit voltage 
considerably once the barrier height is small. In bulk- heteroj unction solar cells based on 
[6,6]-phenyl C^i-butyric acid methyl ester (PCBM) as electron acceptor and poly[2-methoxy- 
5-(3',7'-dimethyloctyloxy)-p-phenylene vinylene] (OC1C10-PPV) as electron donator sand- 



wiched between ITO/PEDO' 
voltage of 0.4 eV is observed 



:PSS and LiF/Al electrodes a reduction of the open circuit 



3Q|. 



B. Steady state 

Now, the steady state situation, where a constant current j flows across the ITO/organic 
semiconductor system, will be considered. The solution is given by Eq. (1251) and the unknown 
constant C% has to be found numerically from Eq.( l28lh 

In Figj2j the solution for the dimensionless electric field F is depicted for current densities 
of j = 10, 100 niA/cm 2 and barrier-free injection (A = 0). 




-8^ — ■ — , — ■ — , — ■ — , — ■ — , — ■ — , 
40 80 120 160 200 



FIG. 2: Distribution of the electric field F in units of Fq = for barrier-free charge- carrier 
injection and a constant current of j = 10 mA/cm 2 (solid line) and j = 100 mA/cm 2 (dashed line) 
as function of the coordinate x in units of ' Itf- 



Comparison with the equilibrium solution shows that the distribution of the electric field 
is, near the electrode, virtually independent of the net current density, but is substantially 
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effected by the current deep in the dielectric bulk. For A = many holes diffuse into the 
organic semiconductor (F' s (0) large) and a strong negative field in the vicinity of the interface 
emerges. This negative field compensates for the diffusion into the organic semiconductor 
to such an extend that the net current reaches j. However, far away from the contact where 
diffusion is negligible, the electric field has to be positive (together with the current). Thus, 
there is a position X = Xq where the electric field changes sign, as can be seen in FigJ21 
Since a vanishing electric field strength is assigned to the ohmic contact itself, the position 
X is often referred to as the virtual electrode [31 ]. 

For a detailed discussion of the field distribution in the insulator, the dependence of the 
position of the virtual electrode, the electric field at the contact F s (0), and the charge- 
carrier density (represented by F^(0)) are shown in FigfS] as functions of A and for two 
typical current densities. Additionally, the distribution of the electric field is depicted in 
FigJUfor different injection barriers and a current density of j = 100 mA/cm 2 . 

For small As, the electric field at the interface is negative, the charge-carrier density 
is large, and the virtual electrode is far inside the insulator. A charge-carrier reservoir is 
formed in the near contact region and the current across the residual insulator X » Xq is 
supplied by this reservoir. From Figj2]andH]it becomes evident that the electric field follows 
F ~ y/X — X for X » X and hence, resembles the field distribution of SCLC assuming 
ohmic boundary conditions 1|. 

Increasing the injection barrier, the position of the virtual electrode is hardly effected 
until the injection barrier exceeds 0.1 eV. For a further increased barrier it approaches the 
physical electrode rapidly. Simultaneously, F s (0) approaches zero and is pinned there for 
quite a wide A-range. However, this does not mean that the ohmic boundary conditions 
(F s (0) = and p s (0) — > oo) assigned to an ideal contact [l| is a good approximation in this 
case. In fact, the density of the charge-carriers at X = is relatively low, though -F s (0) = 



holds. F S (X) reflects neither \/X — Xq nor a constant form. 

By further increasing A, the amount of injected charge carriers remains small, Xq is 
located at the physical electrode and F s (0) becomes positive, while the field in the insulator 
is weakly dependent on X. The charge-carrier reservoir in the insulator is depleted and 
due to the few charge-carriers at the contact, a strong positive field is required to drive the 
current across the near contact region. 

Note that for F s (0) < (or X > 0), the position of the virtual electrode depends 
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FIG. 3: a) Position of the virtual electrode Xq in units of Ij-f for current densities of j = 
10 mA/cm 2 (solid line) and j = 100 mA/cm? (dashed line) as a function of barrier height A. b) 
Electric field strength at X = inside the semiconductor for current densities of j = 10 mA/cm 2 
(solid line) and j = 100 mA/cm 2 (dashed line) as a function of barrier height A. The axis on the 
right hand side indicates the corresponding values for the derivative of the electric field for current 
densities of j = 10 mA/cm 2 (dashed- dotted line) and j = 100 mA/cm 2 (dotted line). 



strongly on the induced current density while F s (0) and i^(0) do not, reflecting the fact 
that close to the contact the equilibrium field distribution is hardly effected by the current 
but further in the insulator it is. For barrier-free injection, the virtual electrode is shifted 
from approximately Xq = 25 nm to 10 nm once the current increases from j = 10 mA/cm 2 
to j = 100 mA/cm 2 . This is illustrated in Fig. [2] and [31 Vice versa, the influence of the 
current density on the field at the contact becomes recognizable for larger barriers, once 
the charge carrier reservoir is depleted and the virtual electrode coincides with the physical 
one. This is due to the fact that a strongly increased positive field is required to support an 
additional current density, since the charge- carrier density is small. 
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FIG. 4: Distribution of the electric field F in units of Fq 



kT 



for a constant current of 



100 mA/cm 2 and barrier height of A = 0.1 eV (solid line), A = 0.2 eV (dashed line), A = 0.3 eV 
(dotted line) and A = 0.35 eV (dash-dotted line) as function of the coordinate X in units oflxF- 

From Figj3]it can be seen that the minimal injection barrier required to result in a match 
of virtual and physical electrode {F s (0) = 0) is shifted to smaller barriers for higher current 
densities. By means of the solution fl25|) . this may be formulated as an exact relation between 
the quantities involved. The requirement F s (0) = results in the particular value of the 
constant C\ = zq(2/\l\) 1 ^ 3 where zq ~ —1.02 is the first zero of the Airy function Ai'(z). 
Using Eq. fT28l) . the current magnitude l , which suppresses the electric field at the interface, 
is determined to 



From the above equation it becomes evident that for increased injection barriers, the minimal 
current density required to obtain a match of virtual and physical electrode is exponentially 
reduced. 

Knowledge about the distribution of the electric field gives access to the voltage drop 
V across the system for a given current density j and hence, to its IV-characteristics. As 
the model consists of infinite semi-spaces, integration over the electric field diverges when 
being carried out over the entire space, i.e. from X = — oo to X = +oo. This holds in 
equilibrium as well as in steady state. However, since the focus of this work is to analyze 
the contact phenomena arising from the conductor /insulator junction, two simplifications 




(37) 
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are introduced. Firstly, according to the definition of the contact potential the voltage drop 
is calculated by introducing a cut-off length L corresponding to the typical thickness of the 
organic layer. Hence, the bulk conduction in an organic layer with finite thickness is taken 
into account. Secondly, since the voltage drop across a conductor bulk of finite thickness is 
small due to its high conductivity, the constant part equal to the asymptotic constant field 
value times the macroscopic conductor thickness is subtracted from the integral over the 
conductor. Then, the voltage drop over the entire system reads: 

K=^l^(0)-^^-21n 

Here, the first two terms in braces represent the voltage drop in the electrode. The third 
term in braces results from the voltage drop within a distance L from the contact inside the 
organic semi-space and the last term corrects the voltage by its equilibrium value, Eq.(|36l). 
The two (not independent) constants F s (0) and C\ result from the steady-state boundary 
condition, Eq. (j28l). 

In Figj5], the resulting IV-characteristics for different barrier heights are presented assum- 
ing L = 100 nm. In the displayed voltage region, the barrier-free contact (A = 0) is able 
to supply more charge-carriers than the bulk of the organic semiconductor can transport. 
Hence, the entire system appears to be space-charge limited with a current approximately 
j ~ V 2 for all voltages. Deviations from the Mott-Gurney law are due to the formation 
of the space-charge region emerging by charge- carrier diffusion at X < Xq. As the current 
increases, the position of the virtual electrode approaches the real electrode and the width of 
this space charge region decreases. Thus, the calculated curve reproduces the Mott-Gurney 
law [l| more exactly for higher voltages. 

In the presence of non-vanishing barriers, the current is substantially reduced. In the 
low voltage regime this reduction is about one order of magnitude once the height of the 
barrier increases from A = to A = 0.4 eV. For high As , one observes a transition from a 
j ~ V 2 to a j ~ V dependence, followed by an exponential current increase with the applied 
voltage. Since the SCLC represents an upper limit for the current flow through the system, 
all IV-curves approach the SCLC regime for even higher voltages. 

The calculated IV-characteristics for high injection barriers can be understood by consid- 
ering some simple approximations. Once the electric field at the interface becomes positive, 
the charge- carrier reservoir established in the organic semiconductor is exhausted and no 
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FIG. 5: IV- characteristics for barrier height of A = 0; 0.2; 0.3; 0.35; 0.4 eV. The curves for A = 
and A = 0.2 eV nearly merge in the above representation. 



virtual electrode is present. Hence, diffusion is negligible all over the space and an electrical 
current is virtually due to charge carrier drift only, 

F' S (X)F S (X) = l. (39) 



The solution depends on F s (0) and reads [32], 



F,(X) = ,JF*(0) + 2lX. (40) 

Analyzing the drift equation (139]) at X = by using the boundary condition, Eq. (j25]) leads 
to an exponential dependence of the dimensionless current i on the local value of the electric 
field 

, = ^(f. F . (0 )-AV (4!) 



a \e c kT J 

where the right hand side of the boundary condition, Eq. (l28]) . has been neglected. Solving 
Eq. (|4T!) numerically with respect to F s (0) for a given current density i determines the electric 
field distribution, Eq. (|40p . and integrating the field over the entire system leads to the IV- 
characteristic. This approximation offers a simple way to calculate IV-characteristics [3^] 
displayed exemplary for a barrier height of A = 0.4 eV in FigEJ 

The agreement with the IV-characteristic calculated using the exact solution, Eq.( [25l) . is 
perfect for all voltages where diffusive transport of charge carriers is negligible. This includes 
the SCLC regime at high bias. Only at low voltages, where the charge carrier reservoir in 
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FIG. 6: IV- characteristic for barrier height of A = 0.4 eV. The solid line has been obtained using 
the exact solution 125\) of the complete drift diffusion equation \19\i. The circles show the IV- 
characteristic calculated with the help of the semi- analytical approach derived for high injection 
barriers, Eqs. lty, and the triangles display the injection current ji n j calculated with the help 

of Eg. (J*fy- The dashed line shows a characteristic current of the virtual electrode appearance. 



the organic semiconductor is not exhausted, deviations are recognizable. According to the 
exact relation, Eq. fl37j) . one can see that below a characteristic current i ~ 4 ■ 10 -3 the 
virtual electrode at X > appears and the diffusion becomes important. 

For a purely injection-limited regime space charge effects are of no importance, resulting 
in a constant electric field all over the semiconductor semi-space. This field coincides with 
the electric field strength at the interface. In such a case F s (0) = VelTF/kTL and Eq.( j4T]) 
governs the injection limited IV-characteristics alone, 

ee s l T FV 



Jinj 



V \r 
eUs—AI exp 

LJ 



A 



e c kTL 



(42) 



34|. Yet, the self-consistent 



Eq. (f42l) resembles the result of the drift-diffusion equation [5|, |33|, 
treatment of the injection problem yields directly a barrier lowering arising from the potential 
energy a charge carrier gains on the conductor side of the interface. In contrast to the 
ln(j') cx W dependence of the Schottky-lowering predicted from the image-charge potential 
in the single-electron picture, the current density depends exponentially on the external 
voltage V. The injection-limited current jinjiY) is depicted in Figj6j 

Once space charge effects become predominant, F s (0) can be neglected in Eq. fj40|) and the 
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temperature-independent Mott-Gurney law is reproduced. This occurs as soon as the field- 
induced barrier lowering has proceeded to such extent that the contact can again establish 
space charge in the organic semiconductor. Hence, a transition from injection-limited to 
space-charge-limited current is observed. As may be seen from FigJHl the crossover from 
injection-limited to SCLC occurs at rather high voltages in presence of a medium injection 
barrier of 0.4 eV. This has its origin in the weak barrier lowering for increasing external 
voltage. Due to the low screening length in the conductor, a high voltage drop across the 
entire system is required to result in an energy gain in the conductor being comparable to 
the injection-barrier energy. Note that in contrast to a space-charge formation resulting 
from diffusion, the electric field at the contact is positive once the space-charge formation 
is solely due to a drift-controlled injection. Such a SCLC regime occurs also for barrier-free 
injection as soon as the diffusive filling of the space-charge region is weaker than its depletion 
due to the external applied field. This, however, requires a very high bias. As a result, the 
SCLC is directly supplied by an efficient charge- carrier drift out of the conductor. 

From Eq. (142l) and Mott-Gurney law one can estimate the upper limit of the barrier 
required to retain SCLC over the whole voltage range: 

A crit = kT + kT\n( ^_L^^_) . (43) 
\27l TFPoo kTj 1 ; 

Assuming the parameters introduced before, we can predict a critical barrier height of 
0.25 eV. This is in slight contradiction to device models in which a Schottky-type bar- 



rier lowering is assumed 24j, |25j, |35|. Here, one would expect the crossover between SCLC 



and injection-limited current to occur at an injection barrier of A crit ~ 0.35 eV. However, 
up to now the experimental data available predict the crossover to be at a barrier height 
between 0.2 eV and 0.35 eV, being consistent with both approaches. 



IV. CONCLUSIONS 



In the calculation of the charge- carrier transport through insulators, a fundamental ques- 
tion about the boundary conditions generally arises when a charge-carrier injecting interface 
has to be involved. Typically, boundary conditions at the interface are chosen, fixing there 
the charge carrier density and/or the electric field. However, at the conductor /insulator 
contact the system is ill-defined, meaning that especially at the interface the charge carrier 



20 



density and the electric field strongly depend on the condition of the system. 

In this paper, a one- dimensional analytical model describing the charge- carrier transport 
across a conductor /insulator junction was presented, where boundary conditions are defined 
far into the conductor and the insulator, respectively. Here, the influence of the two mate- 
rials on each other is negligible so that they can be regarded as independent. Considering 
the Poisson equation and assuming the electric displacement as well as the electrochem- 
ical potential being continuous over the entire system, the electric field distribution and 
the current-voltage characteristic were derived. The model predicts SCLC, injection-limited 
conduction as well as the crossover between the two limits. In most current regimes, the in- 
fluence of the self-consistent treatment is noticeable. For pure injection-limited conduction, 
an injection current similar to the prediction of the drift-diffusion theory was derived. How- 
ever, due to the consistent treatment of the injection problem in one dimension, a barrier 
lowering different to the one predicted from the three-dimensional image-charge potential in 
the single-electron picture emerges. This injection-barrier lowering results from the potential 
energy gain (or loss) of the charge- carriers in the electrode. Also for SCLC, a deviation of the 
Mott-Gurney law is observable as long as the virtual electrode does not match the physical 
one. Only at high bias the self-consistent treatment has no influence on the IV-characteristic 
and the well known Mott-Gurney law is fully reproduced. 
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